Dissipative Quantum Hall Effect in Graphene near the Dirac Point 
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We report on the unusual nature of v = state in the integer quantum Hall effect (QHE) in 
graphene and show that electron transport in this regime is dominated by counter-propagating 
edge states. Such states, intrinsic to massless Dirac quasiparticles, manifest themselves in a large 
longitudinal resistivity p xx ^ h/e 2 , in striking contrast to p xx behavior in the standard QHE. The 
v — state in graphene is also predicted to exhibit pronounced fluctuations in p xy and p xx and 
a smeared zero Hall plateau in a xy , in agreement with experiment. The existence of gapless edge 
states puts stringent constraints on possible theoretical models of the v = state. 



Electronic properties of graphene has attracted signifi- 
cant interest, especially since an anomalous integer quan- 
tum Hall effect (QHE) was found in this material [3, 0] • 
Graphene features QHE plateaus at half-integer values 
of Hall conductivity a xy — (±1/2, ±3/2, ...)4e 2 //i where 
the factor 4 takes into account double valley and dou- 
ble spin degeneracy. The "half-integer" QHE is now well 
understood as arising due to unusual charge carriers in 
graphene, which mimic massless relativistic Dirac par- 
ticles [3} . Recent theoretical efforts have focused on the 
properties of spin- and valley-split QHE at low filling fac- 
tors @, H, @, 0,11, 01 and fractional QHE . Novel states 
with dynamically generated exciton-like gap were conjec- 
tured near the Dirac point TH 12l flU , [l4j | . Experiments 
in ultra-high magnetic fields 15] have so far revealed only 
additional integer plateaus at v = 0, ±1 and ±4, which 
were attributed to valley and spin splitting. 

The most intriguing QHE state is arguably that ob- 
served at v — 0. Being intrinsically particle- hole symmet- 
ric, it has no analog in semiconductor-based QHE sys- 
tems. Interestingly, while it exhibits a step-like feature in 
<7 xy , the experimentally measured longitudinal and Hall 
resistance [15j (p xx and p xy ) display neither a clear quan- 
tized plateau nor a zero-resistance state, the hallmarks 
of the conventional QHE. This unusual behavior was at- 
tributed to sample inhomogeneity [lol ] and remains unex- 
plained. In this Letter, we show that such behavior near 
the Dirac point is in fact intrinsic to Dirac fermions in 
graphene and indicates an opening of a spin gap in the 
energy spectrum [4| . The gap leads to counter-circulating 
edge states carrying opposite spin 0, [B| which result in in- 
teresting and rather bizarre properties of this QHE state. 
In particular, even in the complete absence of bulk con- 
ductivity, this state has a nonzero p xx > h/2e 2 (i.e. the 
QHE state is dissipative) whereby p xy can change its sign 
as a function of density without exhibiting a plateau. 

We start with reviewing the experimental situation 
near v = 0. Our graphene devices were fabricated as 
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FIG. 1: Longitudinal and Hall conductivities a xx and a xy (a) 
calculated from p xx and p xy measured at 4 K and B = 30 T 
(b). The v = plateau in a xy and the double-peak structure 
in a xx arise mostly from strong density dependence of p xx 
peak (green trace shows a xy for another sample). The upper 
inset shows one of our devices. Temperature and magnetic 
field dependence of p xx near v = are shown in the insets 
below. Note the metal- like temperature dependence of p xx . 



described in Ref. Q and fully characterized in fields B up 
to 12T at temperatures T down to 1 K. These measure- 
ments revealed the behavior characteristic of single-layer 
graphene [l[ . Several devices were then investigated in 
B up to 30 T, where, besides the standard half-integer 
QHE sequence, the v = plateau becomes clearly visible 
as an additional step in a xy (FigfTJ). We note, however, 
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FIG. 2: Excitation dispersion in v = graphene QH state for 
a system with and without gapless chiral edge modes, (a) and 
(b) respectively. Case (a) is realized in spin-polarized v = 
state [J], while case (b) occurs when symmetry is incompatible 
with gapless modes, for example, in valley-polarized v = 
state conjectured in Ref. [l5j. In the latter a gap opens at 
branch crossing due to valley mixing at the sample boundary. 



that the step is not completely flat, and there is no clear 
zero-resistance plateau in p xy . Instead, p xy exhibits a 
fluctuating feature away from zero which seems trying 
to develop in a plateau (FigJTJ)). [In some devices p xy 
passed through zero in a smooth way without fluctua- 
tions.] Moreover, p xx does not exhibit a zero-resistance 
state either. Instead, it has a pronounced peak near zero 
v which does not split in any field. The value at the peak 
grows from p xx « h/4e 2 in zero B (7.5 kf2 for the shown 
devices) [1] to p xx > 45kf2 at 30 T (see inset of FigJTb). 

At this point, the absence of both hallmarks of the 
conventional QHE in these experiments can make one 
skeptical about the relation between the observed extra 
step in a xy and an additional QHE plateau. However, 
the described high-field behavior near v = was found 
to be universal (reproducible for different samples, mea- 
surement geometries and magnetic fields above 20 T) . It 
is also in agreement with that reported in Ref. [H| ■ More- 
over, one can generally argue that the QHE at v = can- 
not possibly exhibit the usual hallmarks. Indeed, p xy has 
to pass through zero because of the carrier-type change 
but p xx cannot simultaneously exhibit a zero-resistance 
state because zero in both p xy and p xx would indicate a 
dissipationless (superconducting) state. 

To explain the anomalous behavior of the high-field 
QHE (FigQ]), we note that all microscopic models near 
the Dirac point can be broadly classified in two groups, 
QH metal and QH insulator, as illustrated in FigJ^l 
Transport properties in these two cases are very different. 
The QH insulator (FigJ2t>) is characterized by strongly 
temperature dependent resistivity diverging at low T. 
The metallic T-dependence observed at v — clearly 
rules out this scenario. In the QH metal (FigJ^i), a pair 
of gapless edge excitations (Figj2^,) provides dominant 
contribution to o~ xx , while transport in the bulk is sup- 
pressed by an energy gap. Such dissipative QHE state 
will have a xx ~ e 2 /h 3> & xy , i.e. nominally small Hall 
angle and apparently no QHE. the roles of bulk and edge 
transport here effectively interchange: The longitudinal 
response is due to edge states, while the transverse re- 
sponse is determined mainly by the bulk properties. 

From a general symmetry viewpoint advanced by Fu, 



Kane and Mele [16[ the existence of counter-circulating 
gapless excitations is controled by Z 2 invariants, protect- 
ing the spectrum from gap opening at branch crossing. 
In the spin-polarized QHE state [4j this invariant is given 
by a z . While for other v = QHE states [H EI EE El 
such invariants are not known, any viable theoretical 
model must present a mechanism to generate gapless edge 
states. 

The metallic temperature dependence indicates strong 
dephasing that prevents onset of localization. To ac- 
count for this observation, we suppose that the mean 
free path along the edge is sufficiently large, such that 
local equilibrium in the energy distribution is reached 
in between backscattering events. For that, the rate of 
inelastic processes must exceed the elastic backscatter- 
ing rate: v- inc \ ^> v c \. This situation occurs naturally in 
the Zeeman-split QHE state Q , since backscattering be- 
tween chiral modes carrying opposite spins is controlled 
by spin-orbital coupling which is small in graphene. 

In the dephased regime, the chiral channels are de- 
scribed by local chemical potentials, (pi$(x), whose de- 
viation from equilibrium is related to currents: 

2 2 

/ 2 = ^2, I = h-h, (1) 



where / is the total current on one edge. In the absence 
of backscattering between the channels the currents ii,2 
are conserved. In this case, since the potentials <pi,2 are 
constant along the edge, transport is locally nondissipa- 
tive, similar to the usual QHE [l7j |. 

The origin of longitudinal resistance in this ideal case 
can be traced to the behavior in the contact regions. 
[Note the resemblance of each edge in FigHK with two- 
probe measurement geometry for the standard QHE.] We 
adopt the model of termal reservoirs [ 1 81 ] which assumes 
full mixing of electron spin states within Ohmic contacts 
(see Fig J3}d). With currents ii, I2 flowing into the con- 
tact, and equal currents 1^^ ~ 5(^1 + -^2) flowing out, 

the potential of the probe is Vp ro bc = ^^i™*"*- Crucially, 
using Eq.([T]), there is a potential drop across the contact, 



(2) 



equally for ip\ and tp% . The voltage between two contacts 
positioned at the same edge (see FigfJ^) is equal to V xx = 
\l, which gives a universal resistance value jij. This is 
in contrast with the usual QHE where there is no voltage 
drop between adjacent potential probes [ttI. . 

The longitudinal resistance increases and becomes 
nonuniversal in the presence of backscattering. It can 
be described by transport equations for charge density 

d t ni + d x <pi = j((p 2 - ipi), 

d t n 2 - d x ip 2 = 7(^1 - ^2), ni=Uiipi, (3) 

where 7 _1 is the mean free path for Id backscattering 
between modes 1 and 2, and V1.2 are compressibilities of 
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FIG. 3: (a) Transport in a Hall bar geometry, Eqs.(|3]). The 
edge states denoted by red and blue carry opposite spin. 
Transport through the bulk is indicated by dotted lines. Volt- 
age probes used to measure longitudinal and Hall voltage are 
shown, (b) Voltage probe in a full spin mixing regime [l8l ] 
measures Vprobo = Note finite voltage drop across 

the probe, Eq.©. 



the modes 1 and 2. In a stationary state, Eqs.© have 
an integral I — tpi — if2 which expresses conservation of 

2 ~ 

current / = \l ■ The general solution in the stationary 

current-carrying state is ip\i(x) = vi°2 — jxl. 

For the Hall bar geometry shown in Fig[3^L, taking into 
account the contribution of voltage drop across contacts, 
Eq.©, we find the voltage along the edge V xx — (-fL + 
1)1, where L is the distance between the contacts. In 
the absence of transport through the bulk, if both edges 
carry the same current, the longitudinal resistance is 



Rxx = (jL + 1) 



2e 2 



Pxx = (w/L)R x 



(4) 



with w/L the aspect ratio. From p xx peak value (Fig[T]) 
we estimate jw » 2.5, which gives the backscattering 
mean free path of 0.4 fj,xsx. The metallic T-dependence 
of p xx signals an increase of scattering with T (FigQ}) 
inset). Similarly, p xx growing with B is explained by 
enhancement in scattering due to electron wavefunction 
pushed at high B towards the disordered boundary. 

An important consequence of the Id edge transport 
is the enhancement of fluctuations caused by position 
dependence of the scattering rate j(x). Solving for the 
potentials at the edge, 



I \ (°) 



i I i(y)dy, 



(5) 



we see that the fluctuations in the longitudinal resistance 
scale as a square-root of separation between the contacts: 

SV XX = J f 2 6j(y)dy, 5R X 

J X\ 

where a > 7 _1 is a microscopic parameter which depends 
on the details of spatial correlation of j(x) . Similar effect 
leads to fluctuations of the Hall voltage which has zero 
average value at v — 0. Assuming that the fluctuations 
of the potential at each edge, described by Eq.©, are 
independent, we estimate 5R xy ~ [h/ e 2 ) \J L j a, where L 
is the bar length. 



These fluctuations manifest themselves in noisy fea- 
tures in the transport coefficients near v = 0, arising 
from the dependence of the effective scattering potential 
on electron density. Such features can indeed be seen in 
p xy and p xx around v = (Fig[Tb). As discussed below, 
away from v = bulk transport becomes important and 
short-circuits the edge. This will lead to suppression of 
fluctuations in p xx and p xy away from v = 0, in agree- 
ment with the behavior of the fluctuations in FigQ}}. 

Another source of asymmetry in voltage distribution 
on opposite sides of the Hall bar is the potential drop on 
a contact, Eq. @. This quantity can be nonuniversal for 
imperfect contacts, leading to finite transverse voltage. 
Such an effect can be seen in p xy data in FigQ]near v = 0, 
where Hall effect in a pristine system would vanish. 

To describe transport properties at finite densities 
around v — 0, we must account for transport in the 
bulk. This can be achieved by incorporating in Eq.Q 
the terms describing the edge-to-bulk leakage: 



d x (fii 
-d x <p 2 
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(6) 



where ^1,2 are the up- and down-spin electrochemical 
potentials in the bulk near the boundary. Transport in 
the interior of the bar is described by tensor current-field 
relations with the longitudinal and Hall conductivities 
cii' 2 \ ai^ 2 ^ for each spin component. Combined with 
current continuity, these relations yield the 2d Laplace's 
equation for the quantities V'1,2, with boundary condi- 
tions supplied by current continuity at the boundary: 



4>.V^ + 4>xV^+9(<^-^) = 0, i = 1,2, (7) 

where n is a unit normal vector. [In Eq.([7|) and below 
we use the units of e 2 /h = 1.] To describe dc current, we 
seek a solution of Eqs.© on both edges of the bar with 
linear x dependence ifii(x) — ipf 1 ^ — Sx which satisfies 
boundary conditions |7|). where the functions ijJip, have a 
similar linear dependence. The current is calculated from 
this solution as a sum of the contributions from the bulk 
and both edges. After elementary but somewhat tedious 
algebra we obtain a relation / = 2£/j, where 



\ M 1 ) / W ~( 2 )/ ( 2 ) 

w w Aw ia xy ' I a xx - a xy ' / a xx 

Pxx Pxx £-rs\/(Txx 1 A/ (T xx 



7 2 l + 9 

(8) 

with w the bar width and A = wgj/(2j + g). The quan- 

(1 2) fi 2) 

tities d X y = <J X y ± g/(2j + g) represent the sum of 

the bulk and edge contributions to Hall conductivities, 

and p x X x ,2) are defined as p x l = a xx /(a X y 2 + a x l 2 ). The 

quantity 7, Eq.©, replaces 7 in Eq.Q. At vanishing 
(1 2) 

bulk conductivity, a xx — > 0, we recover 7 = 7. 

The Hall voltage can be calculated from this solution 
as Vh — + V2 — <pv — ( P2 l ), where (pi^i are variables 
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FIG. 4: Density dependence of transport coefficients p^x, p^y 
and = p xy / (p X y 4- p^x ), G X y — pxy/(pxy 4" Pxx ) for a 
Hall bar (Fig(3}, obtained from the edge transport model (|6]) 
augmented with bulk conductivity (Eqs. (|10p .((5)).((9"ll. see text). 
The peak in p xx at v — is due to the edge contribution, 
shunted by the bulk conductivity away from v — 0. Note the 
smooth behavior of p xy near v — 0, a tilted plateau in Gxy, 
and a double-peak structure in Gxx- 



p xy does not show any plateau in the theoretical curve 
(FigHJ), while G xy calculated from p xy and p xx exhibits a 
plateau-like feature. This behavior is in agreement with 
experiment (Fig[T]and Ref.[lB|). 

To conclude, QH transport in graphene at v = is 
due to counter-circulating edge states. In this dissipative 
QHE the roles of the bulk and the edge interchange: the 
edge states dominate in the longitudinal conductance, 
while the bulk conductivity determines the Hall effect. 
This model explains the observed behavior of transport 
coefficients, in particular the peak in p xx and its field and 
temperature dependence, lending strong support to the 
chiral spin-polarized edge picture of the v = state. 

This work is supported by EPSRC (UK), NSF MRSEC 
Program (DMR 02132802), NSF-NIRT DMR-0304019 
(DA, LL), and NSF grant DMR-0517222 (PAL). 



at opposite edges. We obtain Vh = ££, where 



£ = 2w 



riy (a 4- a 



~(2) 
u xy 



A 4- a£! 



^ u XX v XX 



This quantity vanishes at v 



A<Jxx T A<J X x 
n • (!) 

0, since a xy - 



T (2) 
>xy 



(9) 



and 



&xx = Vxx at this point due to particle-hole symmetry. 

In Fig[4] we illustrate the behavior of the longitudinal 
and transverse resistance, calculated from Eqs.©,© as 



PXX = p xy = £7/2 



(10) 



with = 6, gw — 1 (the omitted contact term ([2]) 
is small for these parameters). Conductivities a X x , 

(12) 

Uxy are microscopic quantities, and their detailed de- 
pendence on the filling factor is beyond the scope of 
this paper. Here we model the conductivities a xx by 
gaussians centered at v = ±1, a x x x 2 \v) = e~ A ( v±1 ^ 2 , as 
appropriate for valley-degenerate Landau level, whereby 



a x \ 2 ' is related to a xx ^' by the semicircle relation [lj 



(x^/V^ T 2) 4- (a X x Z >) 2 = 0. In Figdwe used A = 5, 
however we note that none of the qualitative features de- 
pend on the details of the model. 

Fig Q] reproduces many of the key features of the data 
shown in Fig[TJ The large peak in p xx is due to edge 
transport near v = 0. The peak is reduced at finite v 
because the edge is short-circuited by the bulk conduc- 
tivity. The latter corresponds to the double peak struc- 
ture in G xx in FigQJ We note that the part of G xx be- 
tween the peaks exceeds the superposition of two Gaus- 
sians which represent the bulk conductivity in our model. 
This excess in G xx is the signature of the edge contribu- 
tion. The transverse resistance p xy is nonzero due to 
imbalance in a^y' 2 ^ for opposite spin polarizations away 
from the particle-hole symmetry point v = 0. Notably, 
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